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We develop a systematic perturbation theory for the quasiparticle properties of a single impurity 
immersed in a Bose-Einstein condensate. Analytical results are derived for the impurity energy, 
the effective mass, and residue to third order in the impurity-boson scattering length. The energy 
is shown to depend logarithmically on the scattering length to third order, whereas the residue 
and effective mass are given by analytical power series. When the boson-boson scattering length 
equals the boson-impurity scattering length, the energy has the same structure as that of a weakly 
interacting Bose gas, including terms of the Lee-Huang-Yang and fourth order logarithmic form. Our 
results, which cannot be obtained within the canonical Frohlich model of an impurity interacting 
with phonons, provide valuable benchmarks for many-body theories and for experiments. 


The problem of an impurity interacting with a reser¬ 
voir with a continuous set of degrees of freedom plays 
a fundamental role in our understanding of many-body 
quantum systems. Landau and Pekar famously demon¬ 
strated that electrons in a dielectric medium become 
dressed by phonons forming a quasiparticle termed a 
polaron [HE]. Other examples of mobile impurities in¬ 
clude hehum-3 mixed with helium-4 [3| and A particles 
in nuclear matter [3]. Static impurities give rise to the 
Anderson orthogonality catastrophe and the Kondo 
effect l^. With the creation of two-component atomic 
gases characterized by an unrivaled experimental flexi¬ 
bility, the impurity problem can now be studied more 
systematically and from a broader perspective. While 
focus has mostly been on impurities in a Fermi sea (the 
Fermi polaron) (THS), there have been some experiments 
on impurity atoms in a Bose gas pDl - fTS] . With the recent 
identihcation of Feshbach resonances in Bose-Fermi [iBl - 
EH] and Bose-Bose [SI HO] mixtures, the study of impu¬ 
rity physics in a Bose-Einstein condensate (BEC) with a 
tunable interaction is now within reach. 

The impurity problem provides an ideal setting for 
testing many-body theories, and it has yielded funda¬ 
mental insights for the Eermi polaron m- In the case 
of an impurity atom in a BEC—the Bose polaron—most 
studies have either used mean-field theory to study self¬ 
localization [ 22 I - E 5 ] and time dependence for weak cou¬ 
pling [26j , or employed an effective Frohlich model which 
is solved using various many-body techniques P71l32j . 
The Frohlich model, however, ignores interaction terms 
important even for weak coupling, as we shall demon¬ 
strate. The correct microscopic Hamiltonian has been 
used in a field theoretic approach, selectively summing 
ladder diagrams [33], and in a variational approach [34] . 

Since the impurity problem is so useful as a theoretical 
testing ground, it is important to have a quantitatively 
reliable theory, which can serve as a benchmark for other 
many-body theories and for experiments in the weak cou¬ 
pling regime. Here, we provide such an accurate theory 


by developing a systematic perturbation expansion for 
the impurity self-energy to third order in the impurity- 
boson scattering length a. The small parameter of this 
expansion is a/^, with f being the BEC coherence length. 
Also, has to be small for the polaron to be well de¬ 

fined, where ub is the boson-boson scattering length. We 
derive analytical results for the zero temperature quasi¬ 
particle properties of the polaron, showing that the en¬ 
ergy contains a logarithmic term ln(a*/^)a^/^^, where 
a* ^ max(a, ub). When a = ob, the perturbative expres¬ 
sion for the energy has the same form as the celebrated 
result for a weakly interacting Bose gas [351139] . The 
quasiparticle residue and the effective mass are, on the 
other hand, given by analytic power series up to 
We use a T-matrix approach, a technique which has pre¬ 
viously been employed for the Eermi polaron to sum dia¬ 
grams to high order with Monte Carlo methods [4(11442] . 
Our approach uses the correct microscopic Hamiltonian 
rather than the Frohlich model used in Refs. [371432] . as 
the latter is symmetric with respect to the sign of the 
interaction and therefore only contains even powers in 
perturbation theory beyond the mean-field shift. Also, 
the field theory [33] and the variational approach [34] 
miss terms at second order. 

We consider an impurity of mass m immersed in a BEC 
of particles with mass tob. The Hamiltonian is 

+^ X ^B(<?)aI+qaL-qafc'«fc 

k k,k',q 

+ + y X ^(9)4+qaL-q«fc'Cfc: (1) 

k k,k'.q 

where and removes a boson and an impurity, 
respectively, with momentum fc, ef = A:^/2mB and 
e/c = k'^/‘2m are the free dispersions, and V is the sys¬ 
tem volume. The boson-boson Vb (?) and boson-impurity 
interaction V(q) are assumed to be short ranged, and 
they give the usual zero energy scattering matrices 7b = 
47raB/TOB and %, = 2Trajmr respectively, with irir = 


2 


mB'm/{mB + m) the reduced mass, see the Supplemental 
Material [43] . We work in units where h = ks = 1- 
The BEC is assumed to be weakly interacting such that 
it can be described by Bogoliubov theory, i.e. noa^ <C 1 
with no being the condensate density. We define the 
imaginary time Bose Green’s functions as Gii(fe, r) = 
-(TV{afc('r),4(0)}), Gi 2 (fc,T) = -(r.^{a_fc(r),afe(0)}) 
and G 2 i(fc, t) = —aLfc(O)}), where Tr denotes 
time ordering. The Fourier transforms are 


Gii(fc,z) = 


z - Ek z + Ek 


Gl2(fe,z) = 




w ith ul = i ^k/Ek + l)/ 2 , vl = [ik/Ek - l)/ 2 , Ek = 
V^k - ik = ei+Teno, and G 2 i(fe,z) = Gi 2 (fc,z). 

Here ,0 = i2sT is a Bose Matsubara frequency, s is an 
integer, and T is the temperature. 

Perturbation series .—Our aim is to develop a system¬ 
atic perturbation theory in powers of the impurity-boson 
scattering length for the quasiparticle properties of the 
impurity. To this end, we write down all diagrams for 
the impurity self-energy up to third order in the bare in¬ 
teraction V{q). We then formally replace V[q) with Tv 
in each diagram. Diagrams which contain the simulta¬ 
neous forward propagation of an impurity and a boson 
(the pair propagator), such as the second order diagram 
in Fig. [^a) and the three first third order diagrams in 
Fig-[ 2 |;a), can be thought of as coming from the expansion 
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FIG. 1 . Second order diagrams for the self-energy. The 
upper solid black lines are the normal, Gn, and anomalous, 
Gi 2 and G21, boson propagators. The dashed lines are par¬ 
ticles emitted from or absorbed into the BEC, the bottom 
red lines are the impurity propagators, and the external im¬ 
purity propagators are attached to the red dots. The wavy 
lines denote the impurity-boson vacuum scattering matrix T- 
All diagrams come from ladder-type diagrams similar to those 
arising in Eq.[^ The pair propagator Hn appears in (a) and 
(e), 1122 in (b), ni2 in (c), and 1121 in (d) and (f). 


Tv, the first four diagrams in Fig.only contain scatter¬ 
ing of bosons into or out of the BEC, and they can in 
fact be obtained from the Frohlich model if one replaces 
g ^ Tv hy hand [27]. Figures[^e) and (f) contain vertices 
where both the in- and outgoing bosons are outside the 
BEC, and they are not included in the Frohlich model. 
They are, however suppressed by a factor (noa^Y^^, and 
since we only consider terms to lowest order in rioOg, 
these two diagrams will be ignored. Likewise, we do not 
distinguish between ng and n in the following. 

In total, the second order self-energy is 

E2(p) = noTYlTlnip) + 2ni2(p) + n22(p)]. (3) 


Tip) 


Tv 

1 - Tv^ii{p) 


% + tY nii(p) -I— 


( 2 ) 


For T = 0 and p = (0,0), nii(p), ni 2 (p), and 1122 (p) 
may be found analytically [H] to yield 


of the ladder approximation for the impurity-boson scat¬ 
tering matrix in the BEC |44]. Here, nii(p) denotes the 
pair propagator regularized by subtracting the vacuum 
scattering already contained in Tv ■ We use the shorthand 
notation p = (p, w) having analytically continued to real 
energy z ^ oj + z0+. The perturbative expansion ([^ is 
convergent only for small Tv^ii{p), i.e. small a, whereas 
the full frequency dependence of T (p) has to be retained 
when a is large. This approach, which is detailed in the 
Supplemental Material, yields a perturbation series in a 

for the self-energy: E(p) = Si(p) -|- S 2 (p) -b E^Yp) H-, 

where contains diagrams of order a". The first order 
term is the mean-field energy shift Tn with n being the 
total boson density. We now evaluate the next two terms. 

Second order .—We now evaluate the next two terms. 
Second order . The six second order self-energy dia¬ 
grams are shown in Fig. 1 Fig. s a) is given by noT^fHn 
and—together with Fig. 1 e)—comes from expanding the 
T-matrix given by Eq. ([^ to second order in a inside the 
ladder approximation to the self-energy. The three dia¬ 
grams in Figs. Bb) -(d) are given by ng 7 (^ 022 , noTY'n. 21 , 
and ng 7 (,^ni 2 , respectively, with n 2 i = ni 2 represent¬ 
ing the anomalous propagators and II 22 the particle-hole 
propagator |43j . Apart from ladder summations inside 


„ / s 27rnn 

S2(0) = A(a)-—, (4) 

nir ? 

where a = m/ruB is the mass ratio and ^ = 1/y/EimaB 
is the coherence length of the BEC. We have A(l) = 
8 -\/ 2 / 37 r, and A{a) for a general mass ratio is given ana¬ 
lytically in the Supplemental Material |43j . 

Third order. —We now consider the diagrams for S 3 (p). 
They can be divided into three different classes. The first, 
denoted E^aip), is obtained by inserting first order self¬ 
energies Trio for the impurity propagators in the second 
order diagrams depicted in Figs. Ba) -(d). This yields 
four diagrams which are easily evaluated. As we shall 
see, however, these self-energy insertions are canceled by 
a similar first order shift Tno in the impurity energy, 
which must be inserted in the second order diagrams. 

The second class of third order diagrams consists of the 
eight “ladder” diagrams depicted in Fig. [^a). They are 
easily expressed in terms of the two-particle propagators 
Hy. Using the effective propagator G = Gi = G 2 , with 
Gi = G 11 -I-G 12 and G 2 = G 11 -I-G 21 , they can be reduced 
to the two diagrams shown in Fig.j^b). Their sum is 

Ugb = ng7^^[(nii + ni2)^ + (n22 + ni2)^], (5) 
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FIG. 2 . Third order diagrams for the self energy: (a) E36, 
(b) E36 expressed using the propagator G{p), and (c) E3C. 

where we have suppressed the momentum and frequency 
dependence for notational simplicity. 

The third class of diagrams are those where either the 
first or the last two interaction lines are crossed. There 
are 16 such diagrams, but using the propagator G{p) they 
can be reduced to the four terms depicted in Fig. [^c), 
which constitutes a major simplification. Evaluating the 
Matsubara sums and specializing to T = 0 yields for the 
sum of the crossed diagrams 

E3c(p) =2noT^ fdk\ ~ [nn(fc') +ni2(A:')] 

J [z- Ek-ew 

+ ul-ukVk [n,,(fc') + ni 2 (fc')]j, (6) 

with k' = (p + k,z — Ek). 

Equation ^ is, in fact, ultraviolet divergent, the of¬ 
fending terms being UkVk^ii{k') and u^ni 2 (fc'). As 
shown in the Supplemental Material [15] , these two terms 
give rise to a 1/fc behavior of the integrand for Esc for 
large k, which thus appears to be logarithmically diver¬ 
gent. However, this should not cause us too much worry. 
First, the integrand is well behaved at low momenta, 
where a natural lower cutoff is provided by fc ~ 1/^, 
below which the boson dispersion becomes linear. Sec¬ 
ond, the ultraviolet divergence is a consequence of the 
fact that we have assumed a constant scattering matrix 
7b. Retaining the energy dependence of 7b would result 
in an ultraviolet cutoff ^ I/ob. Likewise, replacing Tv 
by the full energy dependent scattering matrix T in the 
diagrams for Ssc gives a cutoff ~ 1/a. We can therefore 
write 

E3(0) = ln(a70 + (7) 

rUr ? 

with a* = max(a, Ob). For the equal mass case we have 
5(1) = 2/3 — -s/S/tt, and the analytic expression for gen¬ 
eral a is given in the Supplemental Material |33|. Note 
that since it does not contribute to 

the self-energy to the order stated in Eq. Q for o > ob. 


It does, however, contribute to the quasiparticle residue 
and effective mass, as we shall see below. 

Quasiparticle energy. —Having evaluated the 28 third 
order diagrams, we can now present a perturbative ex¬ 
pression for the polaron energy given by the solution of 
E{p) = p^/2m -I- E [p, E{p)\. From Eqs. Q and Q, we 
obtain for p = 0 and T — 0 

^ = + Bia)pn{a*/0 (») 

where fl = 2nnQ^/mr is the mean-field polaron energy 
for a = ^. Equation ([^ gives the polaron energy to 
order ln(o*/7o^/7 ^md is one of our main results. We 
see that the small parameter of the perturbation series 
is a/^. The second order term agrees with that obtained 
using the Frohlich Hamiltonian j^. As we can see from 
Fig.j^b) however, the third order logarithmic term comes 
from scattering events where both bosons are excited out 
of the BEG. These are precisely the processes ignored by 
the Frohlich model, which therefore incorrectly predicts a 
vanishing third order term. On the other hand, at fourth 
order in a the Frohlich model has been shown to also 
have a logarithmic contribution m- 

Interestingly, when a = ob , Eq. ([^ has the same struc¬ 
ture as the famous result for the energy of a weakly in¬ 
teracting Bose gas: schematically 5 ^ a[l -I- -b 

na^ln(na^)] jUj. Since it is difficult to measure the 
bulk energy of a Bose gas, the Lee-Huang-Yang (na^)^/^ 
term [55] was measured only recently |2S|, whereas 
the logarithmic correction [37H39] has never been de¬ 
tected. The energy of an impurity atom has, how¬ 
ever, been measured accurately using radio-frequency 
(rf) spectroscopy and our result therefore suggests 
a way to measure beyond mean-field effects in Bose gases 
including logarithmic corrections for the first time. 

Quasiparticle residue and effective mass. —The quasi¬ 
particle residue is given by Z~^ = I — 7;^E. For zero 
momentum, we obtain to third order 

1 — Z ^ = 7tjE2(0,a;)7^„^-b 5c^E3(0,a;)|g 

= 5tjE2(0,a;)|g -b 9tj[E3h(0, w) -b E3c(0,a;)]|g . (9) 

The second line follows from inserting the first order shift 
w = 7/no into the second order self-energy. When ex¬ 
panding in 7/no, this yields a third order term, which 
cancels the third order diagrams Es^; see the Supple¬ 
mental Material [43] . Contrary to the case of the energy, 
the self-energy term Eg^ contributes to the residue since 
and dui'E^c are of the same order. We can evaluate 
dvjZj 2 and dffZ'ib analytically, whereas cI^jEgc has to be 
calculated numerically. We obtain 

Z-^ = l + C{a)—+D{a)^, (10) 

OBt 

where C{a) and D{a) are given in the Supplemental Ma¬ 
terial |35|. For m = mB, we have (7(1) = 2v^/37r and 
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D{1) « 0.64. Equation (101 explicitly shows that the 
polaron is well defined only for a^/aeC ^ 1- 

For an ideal BEC with ob = 0, we have ^ —>■ oo, and 
it follows from Eq. (§ that there are no corrections to 
the mean-field energy up to third order in a. However, 


in this limit Eq. (10) predicts Z = 0 so that there is no 


well-defined quasiparticle, signaling a breakdown of per¬ 
turbation theory. The reason is that the energy of the im¬ 
purity atom is right at the threshold of the particle-hole 
continuum of the BEC, giving rise to a square root en¬ 
ergy dependence of the self-energy and thus zero residue 
as explained in the Supplemental Material |43j . Equiv¬ 
alently, Landau’s critical velocity c = (47raBn)^/^/mB 
above which the polaron decays through momentum re¬ 
laxation, is zero for a noninteracting BEC. 

The effective mass of the quasiparticle is obtained from 
mjra* = Z(1 -f 2mdpiTi). Following steps analogous to 
the calculation of Z, we obtain 


Tn* _ , a?’ 

— — 1 + —- -I- G(q;)— 

m aB4 


( 11 ) 


where F{a) and G(a) are given in the Supplemental Ma¬ 
terial [331 • Ih® equal mass case, we have E(l) = 
16-\/2/457r and G(l) « 0.37. Our result for F{a) matches 
that of Ref. |29| . 

Plots. —In Fig. [^ we plot the zero momentum po¬ 
laron energy, residue, and effective mass, obtained from 
Eqs. (|^, (10), and ([TT]), in the range —0.3 < a/^ < 0.3 
where we expect perturbation theory to be reliable. As 
was discussed above, one should be careful when ob ap¬ 
proaches zero since the quasiparticle residue vanishes in 
this limit. We have chosen ob/^ = 0.1 and depict the re¬ 
sults in the case of equal masses (a = 1), as well as for the 
mass ratios a = 39/87 and a = 87/39 corresponding to 
the experimentally relevant case of a ^®K-®^Rb mixture. 

Consider first the energy. We have i?(0) < 0 [£'(0) > 0] 
for a < 0 [a > 0], corresponding to the attractive (repul¬ 
sive) branch which are both described within our pertur¬ 
bation theory. The second order term gives a significant 
correction to the energy whereas the third order term is 
very small. This is explained in the inset, which shows 
that the third order expansion coefficient B[a) is much 
smaller than the second order coefficient A(a), except for 
a <C 1, so that the third order term is suppressed even 
when a/^ ^ 1. In fact, B{a) —>■ 0 for a very heavy im¬ 
purity with a ^ oo |43] . We note that the polaron can 
form a dimer with a boson for a > 0. This decay process 
is, however, slow in the perturbative regime considered 
here, since the molecule is deeply bound with a binding 
energy —l/2mra^. 

Consider next the quasiparticle residue Z. Here, the 
third order term gives a significant correction, increas¬ 
ing Z for the attractive polaron and decreasing it for 
the repulsive polaron. As we see from the inset, this is 
because the third order coefficient B(a) is larger than 
the second order coefficient C(a). Finally, we see that 



-0.3 -0.2 -0.1 0 0.1 0.2 0.3 

a/^ 


FIG. 3. The polaron energy E, residue Z and effective mass 
m* as given by Eqs. (j^, (10), and (111 for ob/^ = 0.1. We 
show the first order term (dash-dotted lines), second order 
(dashed lines) and third order (solid lines) for equal masses 
(black lines). For a = 39/87 (red lines) and a = 87/39 (blue 
lines) we only show the third order result. In the insets, we 
show the relevant second and third order expansion coeffi¬ 
cients as a function of mass ratio. 


the third order term gives a large contribution to the 
effective mass, decreasing (increasing) it for the attrac¬ 
tive (repulsive) polaron. This is consistent with the inset 
depicting the expansion coefficients F{a) and G(q;). We 
have F{a) —?► 0 and G(a) —?► 0 for a —>• oo [33], indicating 
that the effective mass equals the bare mass for a very 
heavy impurity as expected. 

Varying ^/ub changes the slope of m* and Z as a func¬ 
tion of a/^, but the results are qualitatively the same as 
those depicted in Fig. [^ Varying the mass ratio a also 
changes the relative weight of the second and third order 
terms as explained above. In the Supplemental Material, 
we provide the values of A(a),..., G{a) for a = 0 and 
a —>■ 00 . Intriguingly, G(a), D{a), F{a), and G(a) all 
diverge for a —>■ 0, indicating a breakdown of perturba¬ 
tion theory. In this limit, the atoms in the BEC are much 
heavier than the impurity, and it would be interesting to 
examine how this breakdown is related to the problem of 
a mobile impurity interacting with static scatterers. 

Conclusion. —We developed a systematic perturbation 
theory for the quasiparticle properties of an impurity par¬ 
ticle in a BEC. Analytical results for the energy, residue, 
and effective mass were derived, and the energy was 
shown to contain a third order logarithmic term, whereas 
the residue and the effective mass are given by analytic 
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power series in a up to third order. When ob = a, we ob¬ 
tained the same form for the energy as that of a weakly 
interacting Bose gas, which opens up the possibility of 
detecting corrections to mean-field theory of the Lee- 
Huang-Yang and even a fourth order logarithmic type 
for the first time, using rf spectroscopy. The effects of 
a mass difference between the impurity and the atoms 
in the BEC were analyzed throughout. By deriving rig¬ 
orous results for the quasiparticle properties of the Bose 
polaron, our theory provides a useful benchmark for ap¬ 
proximate many-body theories and for experiments. 
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J. Arlt for useful discussions. R.S.C. and G.M.B. would 
like to acknowledge the support of the Villum Foundation 
via Grant No. VKR023163. 
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Supplemental Materials: Quasiparticle Properties of a Mobile Impurity in a 

Bose-Einstein Condensate 


1 


RENORMALIZATION OF THE CONTACT INTERACTION 


We relate the boson-impurity scattering length to the short-range potential V{p) as follows. The Lippmann- 
Schwinger equation yields the vacuum scattering matrix at vanishing energy in the center-of-mass frame 


% 


1 

y(o)-i-n„’ 


(SI) 


in terms of the Fourier transform of the potential at zero momentum and the vacuum pair propagator n„ = J dfc ^ 2 / 2 m • 

Here we have defined dk = k/(2 tt)^. At the same time, the vacuum scattering matrix is related to the boson-impurity 

scattering length via 


% = 


27ra 

nir 


Similarly for the boson-boson interaction we have 


dTroB 1 

wb 1^(0)“^ — Hb ’ 


(S2) 


(S3) 


with UB = Jdkj^^. 


PERTURBATION THEORY IN THE SCATTERING LENGTH AND THE PAIR PROPAGATOR 


As explained in the manuscript, to obtain a perturbation theory in the scattering length a, we first write down 
the diagrams in increasing order of the bare interaction V{q). We then formally replace V{q) with Tv everywhere, 
which in addition to the obvious substitution V{q) T in the expressions for the diagrams, has one more effect: 
In diagrams where the “bare” forward pair propagator appears when expanding in V{q), we use the regularised pair 
propagator Hu instead when expanding in Tv That is, when making the replacement V{q) ^ Tv, we also make the 
replacement 




dkGii{—k, —ioJv)G{k +p,ioJi, + z) = 
w|(l + /fc) , 


dfc 


nii(p) = dk 


ul{l + fk) 

_Z Ej^ ^k+p 

vlfk , 2771 ,. 


vlfk 


Eh ^k-\-p _ 


Z Eh ^k+p 


Eh Cfc+p 


fc2 


(S4) 


Here G(fc, z) = 1 /( 2 : — e^) is the non-interacting impurity Green’s function, iuju = i{2ii -|- 1)T is a Fermi Matsubara 
frequency, p = {p,z), fh = [exp{Eh/T) — 1]“^ is the boson distribution function. The term 2mr/k^ in Eq. (S4) 


subtracts the contribution already contained within the vacuum scattering, and acts to regularize the pair propagator. 


It comes from inserting the Lippmann-Schwinger equation (SI) in the ladder approximation for T(p) |S44j . 


For T = 0 and p = (0, 0), Eq. (S4) can be calculated analytically yielding 


nii(O) = 


tob 
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1/27^2^ 1 -|- a 


l-^/(a) 

1 -I- a 


(S5) 


This gives nii(O) = uib 12,y/2n'^^ for a = 1. 
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Anomalous and particle-hole propagators 

The anomalous and particle-hole propagators are 

T fk^ 


ni2(p) = -^ X! -iuj„)G{k', iu}^ + z)= Jdk 

OJi, 

n22(p) = -^ X! ^22(-fe, -iu}^)G{k', iuj^ + z) = J dk 


'^k'^kfk 


Ek + efc/ - z Ch' - Ek- z 

uifk ^fc(i + /fc) 

z + Ek-ek' z- Ek- ek' 


where k' = k + p, and 1121 (p) = ni 2 (p). For T — 0 and p = (0,0) we obtain 

TT tr r/ \ 

ni 2 ( 0 ) - . /(a), 

V27r2^ Id- a 

[!-/(«)]• 


(5 6 ) 

(57) 


(S 8 ) 


These expressions are well defined in the equal mass limit a —>■ 1, where we have ni 2 ( 0 ) = m^Bj and 1122 ( 0 ) = 

With Eqs. (S5) and (S 8 ), we have analytical expressions for all pair propagators for p = 0, which is what we need 
to obtain analytical results for the polaron self-energy at zero momentum. 


LOGARITHMIC DIVERGENCE OE Esc AND B{ol) 


Analysing the divergence of £ 3 ^, we investigate the integrand in eq. (§ in the main manuscript for fc —>■ 00 . We 
have Uk -d- 1, Vk ^ TBUome/fc^, Ek -d fc^/2mB, and Ilii{k^—Ek) -d rr?J'^ j^rtiB -f fc^/2M / for k -d 00 , 
where M = m itib- This gives for p = (0,0) 


and 


where we have used 
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(Sll) 


Using these limits in eq. Q, we see that the integrand goes as 1/fc for k ^ 00 . Setting the lower and upper limits of 
the remaining fc-integral to 1 and 1 /a* respectively, we obtain eqs. 0 and (|S13[). 


COEEEICIENTS FOR THE ENERGY: A(q;) AND B{cx) 


We find 


A(a) 


2^2 1 
TT 1 — a 



2a^ 

1 -I- a 



(S12) 


with /(a) = A^/(crZT)7(cr77^arctanwith the definition = i. Note that this expression is 

well defined in the limit a —>■ 1 corresponding to equal masses, where we have A(l) = 8 -\/ 2 / 37 r. 

For the logarithmic term, we find 


77(a) = (1 -I- a) 1 


— arctan \lo? E 2a 

TT 


2-\/a^~+”2a 
7r(l -I- a) 


(S13) 


In the case of equal masses, we have 77(1) = 2/3 — -x/S/tt. 
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COEFFICIENTS FOR THE RESIDUE: C{cx) AND D{a) 


The contributions to the residue from the second order diagrams di^T,2{0, w)|o and the third order “ladder” diagrams 
i9(^S3b(0, w)|o can be calculated analytically, while the contribution from the “crossed” diagrams has to be evaluated 
numerically. We obtain 


with 
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a dimensionless integral, and 


E, = ^p^p^ + 1), = 


Thus, the coefficients for the residue given by Z ^ = 1 + C{a)a‘^/aB^ + D{a)a^ are 
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lUa). (S19) 


COEFFICIENTS FOR THE EFFECTIVE MASS: F{a) AND G{a) 


As for the residue, we can calculate the contributions to the effective mass from the second order diagrams 
dp 2 Jl 2 {p, 0 )[o and the third order “ladder” diagrams dp2Jlob{p,0)\o analytically, while the contribution from the 
“crossed” diagrams has to be evaluated numerically. We have 


2 mdp 2 Y, 2 ip, 0)|o = 
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FIG. SI. Third order diagram for the self-energy obtained by a self-energy insertion in the second order ladder diagram. 
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Thus, the coefhdents for the effective mass given by ^ = 1 -|- F{a)a? JaB^ + G{a)a ^/are 

v/2 3(a + l)-2(a2 + 2)/(a) 

37 r (a -b l)(a — 1)^ ’ 

, 8 (1-b a)^(2-b a^) - 3(a-b l)(3a^-b l)/(a)-b 2a^(a^-b 5)[/(a)]^ 

^“^“3^ (a-l)(a2_l)2 

As stated in the main article the coefficients are well defined in the limit a —>■ 1 and the values are stated in table HI 


EXPANSION COEFFICIENTS FOR a = 0, a = 1 AND a = oo 
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C{a) 

D{a) 

F{a) G{a) 
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C27r 
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TT 

0 0 


TABLE I. Limiting values of the expansion coefficients for the energy, residue and effective mass. 


CANCELLATION OF THE THIRD ORDER SELF-ENERGY INSERTION DIAGRAMS 

To illustrate the cancellation of the third order diagrams obtained by self-energy insertions in the second order 
diagrams, we give here a specific example. Consider the third order diagram depicted in fig. |S1[ which is obtained 
by inserting a self-energy correction in the second order ladder diagram shown in Fig. Ba)- For zero momentum and 
temperature, this gives 
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nlT^ 


dfc 


{uj- Ek- CkY' 


(S25) 


When solving E = E(i?) for the impurity energy to third order in 71,, it is enough to set w = 0 in Eq. (S25), whereas 


we have to insert the first order mean-field energy shift uj = uqTv in the second order diagrams. Inserting this energy 
shift in the diagram depicted in Fig. [^a) and expanding yields 


noTy / dk 


noTi, Ek Ck 


— noTy 




1-f 


Ek + Cfc 


(S26) 


This explicitly demonstrates that the third order term in Eq. (S261 cancels the third order diagram given in Eq. 


(S25). Note that we have ignored the regularising term in the pair propagator in Eq. (S26), since it is irrelevant for 


the present purpose. 


IMPURITY SELF-ENERGY FOR AN IDEAL BEG 

For a non-interacting BEG, there is only one impurity Greens function, Gn, and we furthermore have = 1 and 
Vk = d for the coherence factors. The pair propagator with p = (p, w) becomes the vacuum pair propagator given 
by [SB] 

nvac(p) = (S27) 

\2tt 

and the anomalous and particle-hole propagators both vanish. It follows that the second order self-energy is 

3/2 

S2(p) = -ino7^ ^ x/w - p2/2M. (S28) 

V Stt 

The yTu dependence means that the residue vanishes for a zero momentum/energy impurity, and there is therefore 
no well-defined quasiparticle. 

Furthermore, all third order diagrams vanish forp = (0,0). The ladder diagrams depicted in fig.j^a) are zero as they 
involve the vacuum pair propagator at zero energy/momentum, and/or the anomalous/particle-hole propagators. The 
crossed diagrams in fig. [^c) are also all zero: Either they contain an integral over the particle-hole or the anomalous 
propagator which are both zero, or they contain the pair propagator multiplied by a hole propagator, which is zero 
for a non-interacting BEG. So all 2nd and 3rd order terms vanish for p = (0,0) for a non-interacting BEG, meaning 
that mean-field theory is exact to third order for the energy of a p = (0,0) impurity. 


















